Chaos— induced coherence in two independent food chains 
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presented in Sec. IIIII The type of synchronized regime is 
also investigated in that section. Our main conclusions 
are summarized in Sec. IIVI 
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Coherence evolution of two food web models can be obtained under the stirring effect of chaotic 
advection. Each food web model sustains a three-level trophic system composed of interacting 
predators, consumers and vegetation. These populations compete for a common limiting resource 
in open flows with chaotic advection dynamics. Here we show that two species (the top-predators) 
of different colonies chaotically advected by a jet-like flow can synchronize their evolution even 
without migration interaction. The evolution is charaterized as a phase synchronization. The phase 
differences (determined through the Ifilbert transform) of the variables representing those species 
show a coherent evolution. 
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I. INTRODUCTION 

When dealing with transport processes in complex flu- 
ids flows, the concept of turbulence comes to mind. In 
bidimensional flows, it is possible to obtain a situation 
where chaotic trajectories can be generated by a simple 
and regular velocity field. This situation is called chaotic 
advection 1, 2J. Here we show that chaotic advection in 
a oceanic jet flow can induce coherence evolution in two 
chaotic systems. As an illustrative example, this simple 
and robust mechanism is examined using an ocean food 
chain advected by mesoscale eddies in the ocean. 

In most natural habitats, numerous competing species 
are able to coexist, while in general these communities 
are limited by only few resources (niches). This fact 
contradicts the classical theoretical and empirical stud- 
ies predicting competitive exclusion of all but the most 
perfectly adapted species in relation to each limiting fac- 
tor. Recent developments in the field of chaotic ad- 
vection in hydrodynamical/environmental flows encour- 
age us to revisit the population dynamics of competing 
species in open aquatic systems. A typical model that 
takes into account species interactions is a trophic web 
food chainHHIlll, among them we choose a trophic 
web food chain with a complex behavior. The complex 
behavior, that is, a local disorder, is a requirement be- 
lieved to be necessary for observing non trivial collective 
behavior 0. 

The article is organized as follows. In Sec. we intro- 
duce the ecological model used as well as the flow model 
along with the parameters chosen. The results of the co- 
herent evolution of species inmersed in a chaotic flow are 



II. MODELS 

Simple models for three-species food chains exhibit a 
broad range of non-equilibrium dynamics, from charac- 
teristic natural cycles to more complex chaotic oscilla- 
tions H E S H 113, Two chaotically oscillating 
food web models coupled diffusively may also synchronize 
|l2l H^ . This synchronization phenomenon in coupled 
chaotic systems have been extensively studied .14. .15,] . 
Those systems can display different degrees of synchro- 
nization, namely complete synchronization, phase syn- 
chronization, lag synchronization, and generalized syn- 
chronization |lq . Synchronization by periodic external 
actions in the presence of noise or noise-induced 
has also attracted considerable interest. Recently the ef- 
fect of stirring of chaotic advection in an inhomogeneous 
oscillatory medium was investigated |18|| . In the present 
work we study instead, two chaotic oscillators coupled 
through the chaotic advection of the flow they are im- 
mersed in. 

We use a simple model of two three-species food chains 
immersed in a meandering jet flow. The flow, which is 
laminar and unsteady, produces chaotic advection. One 
important consequence of chaotic advection is the ex- 
ponential separation of initially nearby fluid elements. 
The spatio-temporal dynamics of the two colonies of food 
chains embedded in the time-dependent incompressible 
flow is described by the advection-reaction equations 
which, in a Lagrangian representation, take the form: 
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dr 



V(f,t) 



(1) 
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dt -^^^^-"^ 



f(t)),i,j = l,2 



(2) 



where the second set of equations describes the dy- 
namics of the concentration or the amount of species 
Ui = {ui^Vi,Wi\ contained in a fluid parcel that is be- 
ing advected by the flow described by the first equation 
The flow is assumed to be imposed externaUy, 
so that the population dynamics has no influence on the 
velocity field. At scales large enough (i.e. like ocean cur- 
rents of K, IOO/cto) diffusion effects can be neglected [20^ . 
The coupling between the flow transport capacity and 
the population evolution appears through the spatial de- 
pendence of the Fi{Uj,r) functions. Fi varies from point 
to point, in a fluid element that moves with the flow ve- 
locity. Those functions are evaluated at the position of 
the fluid element at time t, that is, at r(i). 

The population dynamic represented by the function 
Fi(Uj,r) is a metacomunity explicitly modeled by two 
trophic food chains. Standard three level "vertical" food 
chains evolve in every parcel of a well-mixed fluid. The 
resources (i.e. nutrients) ui^2 are consumed by vi^2 (i-G- 
phytoplankton), which in turn are preyed on by top 
predators wi^2 (i-e. zooplankton) . The coupled differ- 
ential equations for the biomass of the different species 
are: 



— = a{ui - UQ{r)) - aifi{ui,vi) 



dui 

~di 

dvi 

~di 
dwi 
~dt 



(3) 



— = -fel-Wl -I- Q;i/i(mi, Wi) - Q!2/2(wi, Wi) , (4) 



— — = -c{wi - W*) + a2f2{vi,Wi) . 



— = a{u2 - uo{r)) ~ aifi{u2,V2) • 



dU2 

~dt 

dV2 

~dt 

dW2 

~dt 



(5) 



(6) 



— = + ai/i('«2,f2) - a2/2(w2, W2) , (7) 



— — = -C(W2 - W*) + a2/2(w2,U'2) 



(8) 



The coefficients a, bi^2, c, represent the respective net 
growth rates of each individual species in the absence of 
interactions among them (ai =0:2 =0). Each 3-species 
model, in the absence of interactions among them, has 
equilibrium or steady state populations (^i 21 "^^i 21 '"'i 2) 
which are the solutions of dui^2/dt = 0, dvi^2/dt = 0, 
dwi^2/dt — (respectively). A linear stability analysis 
yields that the steady state (it* 2 = 0, 2 = 0, wl n_= 0), 
for the chosen parameters (after Blasius et al. [ij) is 
a saddle-node point. We set the origin of each of the 
3-species model as the steady state (in the absence of 
interactions among them and uncoupled with the flow) 
ul 2 — 0, vl 2 — 0, wl 2 = w* > . From the population 
dynamics point of view, this steady state means that the 
predator w is allowed to maintain a low equilibrium level 
even when the prey v, is rare. In other words there are al- 
ternative food sources available for the predator w. The 



two colonies can be distinguished by a parameter mis- 
match of A/ ^ 62 — 61 . The functions fi describe interac- 
tions among the species with strengths a^. We use stan- 
dard interactions of HoUing type II {fi{u,v) = yipf^l) 
describe the competition among species u and v. e in- 
teraction among species v and w is modeled by a Lotka- 
Volterra interactions {f2{v,w) — vw). Equations Q and 
© describe the evolution of (1*1,2) , with net rate a, to- 
wards a space dependent value, ito(r). This term, uo{r), 
is the only explicitly non-homogeneous term, it repre- 
sents a spatially dependent resource (nutrient in a plank- 
ton model) input which could arise naturally from a vari- 
ety of processes such as localized upwelling, river run-off, 
translated as a source or a sink in the flow model. 

The two colonies are chaotically advected by a two- 
dimensional flow. The velocity held of the flow was as- 
sumed to be time dependent, which ensures efficient mix- 
ing. Different flows have proved to produce good stirring 
effect in particles, chemical reactions and plankton (see 
for example and references therein). To illustrate 
this case we choose a flow of geophysical relevance, a jet 
flowing eastward with meanders, of amplitude B{t) and 
wavenumber k in the North-South direction with a phase 
velocity Cx The cartesian components of the flow 

V = {—dtp/dy, dip/dx) are expressed, in nondimensional 
units, in terms of the stream function ^ 



tp{x, y) = 1 — tanh - 



y — B{t) cos k{x ~ Cxt) 



(1 + k^B{tY sin^ k{x - Cxt)) 



1/2 



(9) 



The meander amplitude B{t) is a time-dependent oscil- 
lation, B{t) = Bo + ecos{ujt + 9). 

This flow, representing an open flow, advects eastward 
most of the fluid particles, all together with the species 
contained in each parcel. The source (or sink) of re- 
sources (nutrients) uq{y) is localized at the origin of co- 
ordinates, according to 

u,{x,y)^[^^^ '^^(^^ ^' ^ ^ 

1 elsewhere, 

where the amplitude A is constant. 

III. RESULTS 

The evolution of the colonies in the flow is integrated 
numerically according to the method proposed by Ot- 
tino 2] and later used by others 0, [23 ■ The two- 
dimensional physical space accessible to fluid particles 
is subdivided into regions characterized by different La- 
grangian behaviors. The model we use, without the spa- 
tial dependence, was shown ,13] to have synchronized 
behavior among the top predators species of the two 
colonies when migration, of rate D is allowed. Actu- 
ally in the absence of migration _D = 0, Blasius et al. 
[12. [13 showed that the two colonies would normally be 
nonsynchronized. This nonsynchronizcd behavior can be 
observed in most of the time evolution of the top preda- 
tors, as can be seen in figure ^ In this figure we show 
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FIG. 1: Chaotic time series of top predators, wi and W2 
(dashed line) in the web model. The top images correspond 
to zooms of the main image. The image on the left shows 
the nonsynchronized evolution when the parcel goes through 
a region where the nutrients are homogeneous. On the con- 
trary, when the parcel enters into a region where the nutri- 
ents are non-homogeneous (right) the top predators evolve 
in synchronicity. Parameters for the web food model are 
a = 1,61 = 1.1,62 = 1.055, c = 10,fci = 0.5, ai = 0.2, 02 = 
1.0, w* = 0.006, ui,2(0) = 5.0,^1,2(0) = 5.0,wi,2(0) = 0.0. 
For the flow Bo = 1-2, e = 0.3, uj = 0.4, = 0.12, A = 0.2. 

the temporal evolution of the wi^2 for a fixed parcel in 
the Lagrangian point of view. 

We observe in Fig. ^ that the two web chains are not 
synchronized at the beginning. Then, when the parcels 
enter into the region where the nutrients are spatially 
non-homogeneous, the two colonies start to evolve syn- 
chronously. In the left-top corner a zoom of the first 
time interval is shown. For this period of time the fluid 
parcel goes through a region where the nutrients are ho- 
mogeneous. A nonsynchronized evolution is observed, as 
was expected for these parameter values |T^ . A zoom 
of the synchronous evolution can be seen in the right- 
top corner. This inset corresponds to the interval when 
the spatial dependence, that is, the non-homogeneous 
sources or sinks, advects the resources (nutrients). This 
figure shows also that the evolution of the two subsys- 
tems undergoes a transition to another chaotic attractor, 
adjusting their rhythm due to the interaction of the flow. 

The coherent evolution of the top predators of the 
two colonies wi and W2 can be explicitly shown plot- 
ting 'W2 vs .wi as in figure El In the top panel (Fig. El a) 
) W2 vs. wi is plotted during the time where the source 
is not forcing the system. The clouds of points clearly 
shows the uncorrelated behavior of the two variables. In 
Fig. Elb), on the contrary, W2 and wi display a coherent 
evolution. 

It is well-known that two chaotic systems could dis- 
play different synchronized regimes [Illlllll. Thus, 
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FIG. 2: Projections of the phase portrait on (uii,ui2) plane. 
Panel (a) and (b) correspond to the cases where the nutrients 
are homogeneous and inhomogeneous respectively. Notice the 
perfect synchronization. Parameters for the web food model 
and the flow are the same as for Fig0 



we investigate the regime of synchronization of the two 
colonies, as well as the influence of the degree of mixing 
power of the flow. 

A closer look at the synchronized evolution reveals that 
this is not a complete synchronization regime. In fact, 
wi ^ W2 ~ bi — &2 (where — means a temporal average 
of ...). Another possible scenario is that the two colonies 
are in a phase synchronization regime. 

To describe the phase synchronization, we need to in- 
troduce corresponding quantities. The phase of the sig- 
nal (the time evolution of the population density of one 
of the species) can be obtained by different ways. We 
calculate the phase, using the standard construction of 
the analytic signal [isl l2l |. Complex signals are ob- 
tained from the real signals are rewritten as a complex 
signals zi,2 = wi,2(0 + *^[w'i,2(i)] = |ai,2(i)|e*-^(*), be- 
ing H[wi^2{t)] the Hilbert transform of wi,2(i) ■ 

In figure we plot the relative phase difference 
Auj{t) = 01 (t) — (/'2(t) as a function of time. The method 
of phase estimation has several advantages and some 
drawbacks (see for a more detailed discussion and 
Cap. 6 and Appendix 2 of ^3)- The evolution of the 
point in the complex {wi^2(t), H[wi_2it)])-pl£L'ne rotates 
around two different centers of the two chaotic attrac- 
tors. The density of population wi^2 evolve in a chaotic 
attractor while nonsynchronized, then after a transient 
time they fell into another chaotic attractor (with differ- 
ent signal mean value) and they synchronized in phase. 
After another transient time (where the mean value is 



4 




FIG. 3: Plot of phase the difference Aaj(t) = (jjii't) — <^2(t) as 
a function of time. The regions where the phase difference 
growths (meaning no synchronization) correspond to non- 
chaotic advection while the regions of constant phase differ- 
ence are associated with synchronized evolution. Parameters 
for the web food model and the flow are the same as in Fig. 

m 



neither of the previous one) they rotate once again non- 
synchronously in the first chaotic attractor. The transi- 
tions between the different attractors and the time spent 
in each attractor is clearly reflected in Figure 13 The 
first region shows the two colonies not synchronized, the 
phase difference Auj{t) grows with time. The next win- 
dow time corresponds to the transient time when the 
colonies abandon the nonsynchronized chaotic attractor, 
the phase can not be estimated with the above procedure. 
Once the parcel reaches the region where the source (or 
sink) of nutrients is located, (uo(r) ^ 0), the phase differ- 
ence between the two patches drops to a constant. The 
two species are synchronized. The two species evolve in 
another chaotic attractor with well defined mean value. 
Then, as the parcel leaves the upwelling region (or sink), 
there is a transient during which the signal leaves the at- 
tractor evolving to the nonsynchronized one. Once again, 
the phase difference in the transient time the phase can 
not be calculated. After that transient the evolution of 
the two species is once again nonsynchronized. 

The flow can influence the coherence evolution of the 
two colonies. Z. Neufeld and coworkers 03| has shown 
that oscillators advected chaotically by a flow can pro- 
duce collective oscillations or oscillator death by control- 
ling the mixing capacity of the flow, actually the stirring 
rate of flow. It has been shown ^2T\ that the mixing ca- 
pacity of the type of flow we are using in this work can 
be modified by three parameters that govern the time- 
dependent oscillation of the meander amplitude, namely 
Bq, e, uj. We choose the values used by Cencini et al. 22] 
and later by Lopez and coworkers jlOj originally moti- 
vated mainly by observations in oceans jets. These are 
the critical values for obtaining "large scale chaos". It 
is under this situation that exchange of particle between 
north-south part is more favorable (more mixing). The 



parameters Bq = 1.2, e = 0.3,0; = 0.4 were chosen by 
Cencini et al. to be greater than the critical value in or- 
der to have great power of mixing. Different collective 
behaviors of the two colonies are expected when the mix- 
ing capacity is changed through any of the three param- 
eters involved, as they will be changing the mean value, 
the amplitude and the stirring rate of the inhomogene- 
ity. Different situations may arise though, changing the 
parameters of the flow (as well as the parameter of the 
colonies dynamic), these situations are discussed some- 
where else 25]. 
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FIG. 4: The two top predators wi^2 and the flow forcing uo(r) 
as functions of time are shown. It is shown only the region 
where the nutrients are non-homogeneous. The two panels 
show the independent behavior of the variables wi,2 of the 
dynamical system related to the forcing {lu — 0.5,^ = 0.5). 
Parameters for the web food model are the same as for FiglTI 

We would like to remark that the two population 
{wi,W2) are in fact synchronized. That is they are not 
just two (chaotic) oscillators passively following the same 
forcing. It can be observed in Fig. 0] that the top- 
predators zi'i_2 are not adiabatically following of the forc- 
ing. Indeed for the values where the flow keeps constant 
(and different from zero) the variables wi_2 are oscillat- 
ing. 



IV. CONCLUSIONS 

In summary, we have addressed the evolution of two 
food web models immersed in a flow. The chaotic advec- 
tion, a mixing mechanism present in the real ocean, was 
shown to induce a coherence evolution of two species of 
different colonies. In particular we have considered two 
food web models of a three-species food chain, each ad- 
vected by a jet-like flow. The population model consid- 
ered here represents a quite general population dynamic, 
and the main issues found here, namely the possibility 
of finding coherence evolution of two species, may be 
present in biological transport situations. There still re- 
main open questions in this issue, such as mapping in the 
phase parameter space all the possible collective behavior 
of the two colonies and the transition between the differ- 
ent attractors. A simpler dynamical model (although less 
ecologically plausible) may help to examine the influence 
of different flows. Finally, we stress that the coherence 
evolution of two species as a result of a mixing property 
of the flow they are immersed in is a powerful process. 
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It has the potential to shape the distribution and abun- 
dance of aquatic species in a current flow with important 
implications for ecological dynamics in fluid flows. 
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